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Using the finite temperature quasi-particle random phase approximation (FTQRPA) on the basis 
of finite temperature Skyrme-Hartree-Fock + BCS method, we study -decay half-lives for even- 
even neutron magic nuclei with N=82 in a finite temperature environment. We find that the 
decay half-life first decreases as the temperature increases for all the nuclei we study, although 
the thermal effect is found to be small at temperatures relevant to r-process nucleosynthesis. Our 
calculations indicate that the half-life begins to increase at high temperatures for open shell nuclei. 
We discuss this behavior in connection to the pairing phase transition. 



PACS numbers: 21.60.Jz, 25.30.Pt, 26.20.Np, 26.30.Jk 



I. INTRODUCTION 

-decay of neutron-rich nuclei is one of the impor- 
tant subjects for r-process nucleosynthesis. In the r- 
process, nuclei rapidly capture neutrons and reach the 
neutron-rich region, until the timescale of neutron cap- 
ture is comparable to that of the photodisintegration in 
the vicinity of the neutron shell gaps A/'=50, 82, and 126. 
The -decay becomes important mainly at this point. 
It increases the atomic number of neutron-rich nuclei, 
and eventually enables them to go on capturing neutrons 
again toward heavier nuclei. Therefore, the -decay 
half-lives of neutron-rich nuclei determine the r-process 
time scale, and thus considerably influence the final abun- 
dance of elements. Likewise, the /3+-decay plays a deci- 
sive role in the evolution of rp-processes elements [1]. 

Most -decay rates of neutron-rich nuclei relevant to 
the r-process have not yet been measured experimentally. 
Therefore, r-process calculations have to rely on a the- 
oretical estimate of /3-decay half-lives. Several theoreti- 
cal approaches have been developed so far. One of the 
most widely used theoretical methods is the gross the- 
ory [il, which describes the /3-decay rates with a sum 
rule approach supplemented by a statistical treatment 
for final states. Although it has enjoyed considerable 
success, it is not clear how well the shell and pairing ef- 
fects for weakly bound systems are treated in the theory. 
Another approach is the shell model, which successfully 
reproduces the experimental half-lives of waiting-point 
nuclei at 7V=50, 82, and 126 |3,[i,[5|. However, a large- 
scale shell model calculation for a systematical study for 
heavy nuclei along the r-process path has been limited so 
far. 

The proton-neutron quasi-particle RPA (pnQRPA) 
i, 0, H i, H [m, E 111, Q E H [13 is suitable 
for bridging the gap between the two approaches. Being 
the microscopic approach, the pnQRPA properly takes 
into account the shell and pairing effects, and more- 
over it is ideal for a systematic study. The strength 
that contributes to the [3~ decay mainly comes from a 
small low-energy tail of the Gamow- Teller (GT) distri- 



bution, which is in general difficult to reproduce accu- 
rately with the pnQRPA. However, the pnQRPA ap- 
proach based on the microscopic self-consistent mean- 
field framework has successfully reproduced the /3-decay 
half-lives for neutron-rich isotopes by appropriately ad- 
justing the proton-neutron pairing strength in the isospin 
T=0 channel [i, [HI, [13 • 

The r-process takes place in an environment of high 
temperatures (T ~ 10^ K) and high neutron densities 
{p ^ 10^^ neutrons/cm^). In this environment, a part 
of excited states is thermally populated, and in princi- 
ple one needs a finite temperature treatment for /3-decay 
calculations for r-process. Notice that the thermal effects 
affect especially low-lying states, which are important for 
the /3-decay. The thermal effects on the /3-decay rates has 
been studied with an independent particle model [H| 
and with the finite range droplet model (FRDM) plus 
gross theory [19]. The temperature dependence of elec- 
tron capture rates has also been studied with large-scale 
shell model calculations |[5j| as well as with pnRPA with 
Skyrme interaction(2Ql|. 



In this paper, we assess the thermal effects on the /3- 
decay of neutron-rich nuclei using the pnQRP A ap proach. 
A similar attempt has been done in Refs. pisl. fl4 121I. [22| . 
but they have used a schematic separable force for the 
particle-hole interaction. Some of them have neglected 
also the proton-neutron pairing correlation. We instead 
carry out our calculations based on the finite tempera- 
ture Skyrme Hartree-Fock + BCS method, together with 
a contact force for the proton-neutron particle-particle 
interaction in pnQRPA. 



The paper is organized as follows. In Sec. II, we sum- 
marizes the theoretical method for finite-temperature 
QRPA. In Sec. Ill, we show the calculated results for 
the isotones with neutron magic number TV = 82, which 
are relevant to the r-process nucleosynthesis. In Sec. IV, 
we give a summary of the paper. 
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II. THEORETICAL METHODS 

Finite temperature Hartree-Fock + BCS 
method 



with /i = A£^/10. The variable AE is determined so as 
to satisfy, 



(7) 



In order to study /3-decays at finite temperatures, we 
first construct the basis states using the finite temper- 
ature Hartree-Fock+BCS method til, (H]. The for- 
mahsm of the finite temperature Hartree-Fock + BCS 
method is almost the same as that at zero-temperature 
m, [23, HSi] , except for the particle number and pair- 
ing densities. At zero temperature, the single-particle oc- 
cupation probability is given by the BCS occupancy 
vf. On the other hand, at finite temperatures (3 = l/kT, 
k being the Boltzmann constant, it is modified to. 



ni{T)=MT) + tc.nh(^^]vf, 
MT) = (alai) 



(1) 



1 + expiPEi) ' 



where and fi{T) are the creation operator and the 
occupati on probabihty for a quasi-particle, respectively. 
Ei = \J (ci — A)2 + (Aj)2 is the quasi-particle energy, 
where and A are the single-particle energy and Fermi 
energy, respectively, and the pairing gap obeys the 
gap equation, 



"2^ 

j>0 



^7 1 

tanh 



^ 2 



(2) 



Here, V is the pairing interaction and i is the time- 
reversed state of i. 

With the densities obtained with the single-particle oc- 
cupation probabilities n^, the self-consistent solution is 
sought by minimizing the free energy. 



F = E -TS{T), 



(3) 



where E is the Hartree-Fock energy and and S{T) is the 
entropy defined as, 

S{T) = -kJ2MT) In MT) + (1 - /i(T)) In (l - /^(T)) . 

i 

(4) 

In the calculations shown below, we use the smooth 
cutoff scheme for the pairing active space, following Ref. 
(26I . [2§|. That is, the quasi-particle energy and the gap 
equation are modified to Ei = (e^ — A)^ + (7iAi)^ and 



A, = -iyF,,,^^tanh^^ 

j>0 ^ 



(5) 



respectively. Here, the cutoff function is defined as [26|, 



1 



1 + exp [(ei - A - AE)/ij] ' 



(6) 



where Nq is the number of particle for proton (q=p) or 
neutron (q=n). 

In our calculation, we employ the zero-range density 
dependent force. 



V^pair(r - r') = V^^ (^1 



Po 



5{r-r'), (8) 



for the like-particle (proton-proton and the neutron- 
neutron) pairing interactions. We neglect the proton- 
neutron pairing for the BCS calculation, although it is 
taken into account in the QRPA calculation, because we 
are interested in neutron-rich nuclei, rather than N::^Z 
nuclei, in which the proton-neutron pairing plays a mi- 
nor role. For the parameters for the pairing interaction 
in Eq. (|8]), we fix po to be 0.16 fm~^ and adjust the 
strength parameter so as to reproduce the empirical 
values for the pairing gap obtained from the three-point 
mass difference A^^\N + 1) [30]. 



B. Finite temperature quasi-particle random phase 
approximation (FT-QRPA) 

Collective motions of hot stable nuclei have been stud- 
ied with the finite temperature random phase approxima- 
tion (FTRPA) [si, m, [Hi. it was firstly developed 
for studying the giant dipole resonance of a hot com- 
pound nucleus formed in heavy-ion reactions. To discuss 
the property of hot exotic nuclei, the finite temperature 
quasi-particle RPA (FT-QRPA) was recently employed in 
Ref. [35|. The finite temperature proton-neutron QRPA 
has also been developed in Refs. [2l|, [22| for a separable 
interaction. 

The applicability of the finite-temperature RPA has 
been assessed with the Lipkin-Meshkov-Glick method 
m, 113, m, llil. These studies have shown that the 
finite-temperature RPA works satisfactorily well for the 
total strength. When the interaction is small so that the 
ground state is spherical, the FTRPA also yields a rea- 
sonable strength function. Because the proton-neutron 
coupling is usually weak (that is, the isovector interac- 
tion is not large enough to 'deform' the ground state in 
the isospin space), we argue that the finite-temperature 
pnQRPA provides a reasonable tool to discuss the ther- 
mal effects on the /3-decay rate. 

At finite temperatures, the quasi-particle states are 
thermally occupied according to the quasi-particle occu- 
pancy fi{T) in Eq. ([T]). Therefore, the excitations involve 
both two-quasiparticle excitations and one-quasiparticle 
one-quasihole excitations, as is schematically shown in 
Figure [H This can be understood as follows (see Fig. [2j). 
At zero-temperature, the excited states corresponds to 



3 



Ei 



T=0 



Ej 




quasi -hole 



FIG. 1: (Color online) A schematic illustration for the ex- 
cited states in QRPA. At zero temperature (the left panel), 
excited states correspond to two quasi-particle (2qp) states, 
while one-quasi-particle one-quasi-hole (Iqp-lqh) states are 
also involved at finite temperatures (the right panel). The 
occupation probability of one quasi-particle states is given by 
MT) in Eq. ([l]). 




2qp state (k, i) 

Ek-Ej 

qp-qh transition 
2qp state (i, j) 



ground state 



FIG. 2: (Color online) A schematic illustration for excita- 
tion scheme in finite-temperature QRPA. Transitions from 
the ground state corresponds to a two-quasiparticle excita- 
tion, a^a\ whereas the transitions among two-quasiparticle 
states are described by the operator a^a. 



two-quasi-particle states built on the quasi-particle vac- 
uum. At finite temperatures, these excited states are 
thermally populated. The transitions among the two- 
quasi-particle states are then described by the operator 
a^a, e.g., 

ala\\0) = alaj[a\a\m. (9) 
The energy change for this transition is 

= {Eu + Ei) - {Ei + Ej) =Ek- Ej. (10) 

The transition operator at a finite temperature thus 
reads [33], 

a,f3 

(11) 

where a and P run over proton and neutron levels, re- 
spectively. The first and third terms in Eq. ([TT]) cor- 
respond to the transitions among the two-quasiparticle 



states, which vanish at zero-temperature. The QRPA 
equation can be derived from the equation of motion, 
{\[SQ,[H,Q^]]\) = ^qrpa(|[^Q,Q1|), where ^qrpa is 
the QRPA excitation energy and SQ is any one-body op- 
erator. This yields, 



( C a D 

A B 

-D -b -C -a 
\ -b^ -B -a^ -A J 



b\ fP\ 

X 

Q 



£^QRPA 



X 

Q 



(12) 



where the elements of the matrices A^ B^C, D^d, and b 
are given by [33], 



Aa(3a'(3' — \/l - fa - ff3A'^p^,p,y^l - fa' - f(3 

+ [Ea + E^)5 
Ba(3a'(3' = \/^ — fa — f (3^ a(3a' (3' \/ ^ — fa 



Ca[3a'(3' — \/ f (3 ~ f aC a(3a' (3' \/ f (3' — fa' 

+ — E 13)5 oca' ^[3 [3' 

Daf3a'f3' = \/ f (3 — f a^ ai3a' (3' \/ f i3' " fa' 

Q^a(3a'(3 

^a(3a'(3 

with 



(13) 



A' 

^af3a'f3' 
Baf3a'f3' 



^Lf3a'f3' 



Da(3a'(3' 
C^a(3a'(3' 
baf3a'f3' 



— V fp ~ fcxCLa(3a'(3' \/^- fa' — ^/3' 
= \/ f(3 — faba(3a'(3'\/^ " fa' " f(3'-, 

Va(3a'(3'{UaU(3Ua'U(3' + VaVpVa'Vp') 
+ ^aP'pa'{^aV(3Ua'V(3' + VaUpVa' Up^) 
^af3'a'p{'^a^f3^a''^(3' + VaUf^U^'V (3') 

- ^a(3a'p'{'^aU(3Va'V(3' + VaVpUa'Up^) 
Va(3'(3a'{u 

- ya^^'ai^aVf3Ua'V(3' + VaUpVa' Uf3') 

- yaf3a'p'i^aV(3Va'Uf3' + VaUpUa'Vp^) 
~^ '^a^'a' fdi'^c^ U(3Va'V(3' ^VaVpUa'Up') 

^a^a'f3'i^aU(3Va'V(3' - UaVpUa' Uf3') 
~ ^a^'/3a' i^a V(3Va'Ui3' - UaUpUa'Vp') 
Vaf3a'f3'{VaU(3Ua'Uf3' - U^V [3Va' V f^') 

- Vaf3'f3a'{VaVf3Ua'Vf3' - UaUpVa' Up') . 



(14) 



Using the solution of the QRPA equation, the strength 
function S^{E) for the GT transition is calculated as , 



S^{E,) 



1 



1 - exp{-pE^) 



UaU(3Pa(3 + VaV(3Qa(3)Vf(3-fa 



a>(3 



^{UaVpX^^ + VaU(3Y^^)^yl- ff3- fa^ 

xS{Ea-Ep-E,), 



(15) 



where Oqj. = ar^. For T = 0, one can see that the Eqs. 



(|T2|) and ([15]) are correctly reduced to the usual QRPA 
equation at zero temperature. 
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In our calculations, we use the to and ts terms in the 
Skyrme force [40] as the residual interaction, 



v{r^ r') — 



^ + {a-a){T-T)5{r-r% (16) 



for the Gamow- Teller transition. For the particle-particle 
matrix elements (the proton-neutron isospin T=0 pair- 
ing) in Eq. ([T4|), we use a 5-type interaction. 



T4,„(r,r') = 0(r-r'). 



(17) 



We can regard VSn ^ ^^^^ parameter as has been dis- 
cussed in Ref. because we do not take into account 
the T=0 pairing in the Hartree-Fock calculation. The 
Gamow- Teller low-lying strengths are sensitive to the 
T=0 pairing, and we adjust the value of V^^ to reproduce 
the known experimental (5 half-life at zero temperature 

We solve the QRPA equation by diagonalizing the 
QRPA matrix in Eq. (p!2]) . In order to include continuum 
states, we discretize them with a box boundary condition 
with the box size of 15 fm. We include the single-particle 
states up to e = 20 MeV, and truncate the QRPA model 
space at the two quasi-particle energy of E2qp = 70 MeV. 
Our pnQRPA calculation is not fully self-consistent, since 
we do not include all the residual interaction terms in the 
Skyrme functional. We thus scale the residual interac- 
tion Eq. (p!6|) so as to reproduce the spurious translational 
mode (that is, the isoscalar dipole mode) at zero energy 
at every temperature we consider. 



III. RESULTS 



Temperature dependence of GT strengths for 
N=82 Nuclei 



Let us now numerically solve the pnQRPA equations 
and discuss the temperature dependence of the GT 
strengths for even-even N=S2 nuclei, ^^^Sr, ^^^Zr, ^^^Mo, 
^^^Ru, ^^^Pd, and ^^^Cd. For this purpose, we mainly use 
the SLy5 force [41] for the Skyrme parameter set. We 
set the proton pairing strength Vp = —1300 MeV-fm~^ 
so as to reproduce the empirical pairing gap of ^'^^Cd, 
that is, A^p\z + 1) = 0.92 MeV. The proton-neutron 
pairing strength in Eq. (p!7|) is adjusted to V^^ = —360 
MeV-fm"^ so as to reproduce the experimental /3-decay 
half-life of ^^^Cd (0.195 sec.) [42], and use the same value 
for all the other nuclei. 

We find that the strength function is almost the same 
as that at T = for temperatures less than T = 0.2 MeV, 
which is considered to be the standard r-process temper- 
ature at the initial condition |43|. Figure [3] shows the GT 
strengths at T = 0.0 (the sohd line) and T = 0.3 MeV 
(the dashed line) for the ^^^Zr, ^^^Ru, and ^^°Cd nuclei 
as a function of — Ei, where Ei and E^ are the en- 
ergy of the initial and final states, respectively (See FigHl 
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FIG. 3: (Color online) The Gamow- Teller strength functions 
for the ^22zr, ^^^Ru, and ^^°Cd nuclei at T = 0.0 (the solid 
line) and T — 0.3 MeV (the dashed line). These are plotted as 
a function of — Ei, where Ei and E^ are the energy of the 
initial and final states, respectively. The peaks which vanish 
at T = 0.0 MeV are indicated by the arrows. For ^^^Ru and 
^^^Zr, the strengths at T = 0.3 MeV are scaled by a factor of 
16. 



and Eq. ([2T]) ). Those strength functions are smeared 
with the Lorentzian function with the width of 0.1 MeV. 
The strengths at T = 0.3 MeV for ^^^Ru and ^^^Zr are 
multiplied by a factor of 16 for the presentation purpose. 
One sees that some new peaks, indicated by the arrows, 
appear at T = 0.3 MeV, which originate from the tran- 
sition from the excited states. Their strengths are of the 
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FIG. 4: A schematic illustration for the /3~ -decay scheme 
at finite temperatures. The transitions at zero temperature 
are indicated by the thick solid arrows, while the additional 
transitions at finite temperatures by the thin solid arrows. 



order of 0.1 on average, which are approximately of 0.1 
% of the sum rule. Despite its small value, these con- 
tributions to the /3-decay half-life cannot be neglected as 
we will discuss in the next section. 



B. /3-decay Half-Lives 

We next calculate the /3~ -decay half-lives. Since the 
contribution of the GT transition to the total -decay 
rate is much larger than the Fermi transition ^4^], we take 
into account only the former. The /3-decay half-life 
can be calculated with the Fermi Golden rule as 



1 



/2 



In 2 
ln2 h 



(18) 



where is the /?-decay rate. Gp = 1.1658 x 
10~^^ MeV~^ is the Fermi constant, and qa = Ga/Gv 
is the ratio of the vector and pseudo vector constants, 
which we set 1.26. The function p{E,Ee) is the phase 
space factor for the outgoing electron and anti-neutrino 
given by. 



p{E,E,) = ^^XKI^^E - E,fF{Z,E,), (19) 

where Ee is the energy of the electron and Z is the atomic 
number of the parent nucleus. F{Z^Ee) is the Coulomb 
correction factor given by [46], 



r(7 + iv) 



r(27 + i) 



(20) 



where 7 = (1 — Z'^a^Y^'^ and v = {ZaEe/pec), a being 
the fine structure constant, ke = Pe/h is the electron 



wave number and T{x) is the gamma function. The en- 
ergy Ei - E^ in Eq. is related to the pnQRPA 
energy ^^qrpa as 



E, 



AM., 



-H 



RPA 



(21) 



where AMn-H = 0.78227 MeV is the mass difference 
between a neutron and a hydrogen atom. 

Figure [5] shows the /3-decay half-lives normalized to 
that at zero temperature, as a function of tem- 

perature T. In order to check the parameter set depen- 
dence of the Skyrme functional, the figure also shows the 
results with the SkM* parameter set [47]. One sees that, 
as the temperature increases, the /3-decay half-life first 
decreases gradually for all the nuclei we study for both 
the parameter sets. 

One can also see that the temperature dependence is 
the stronger for the larger atomic number. For instance, 
at T = 0.8 MeV, the ratio is around 0.2 for 

^^^Cd both for the SkM* and SLy5, while it is about 0.9 
for ^^^Sr. This behavior is related with the number of 
the GT peaks. Figure [3] indicates that the number of GT 
peaks decreases gradually with the atomic number. This 
is due to the difference between the proton and neutron 
Fermi surfaces. For ^^^Cd, the number of GT peaks is 
only two at T=0, and the thermal effects are relatively 
large. On the other hand, the effects are less significant 
for ^^^Zr because there are already many strengths at 
T = 0.0 MeV. 

For ^^^Zr, ^^'^Mo, and ^^^Ru, the half-lives begin to in- 
crease at temperature around T = 0.6 — 0.7 MeV. This is 
related to the temperature dependence of £^qrpa, which 
also influences the phase space factor in Eq. ([19]). That 
is, when £^qrpa is large, the phase space factor is also 
large, resulting in a large /3-decay rate {i.e., a short half- 
life). Since the excitation energy from the ground state, 
Ef^ = Ei — Eq and ^qrpa now depend on the temper- 
ature, the /3-decay half-life may not behave in a simple 
way as a function of T. 

A simple estimate of the thermal effect on the /3-decay 
half-life can be made by disregarding the temperature 
dependence of S~{Em) in Eq. ([18]) except for the the 
thermal population probability of excited levels. We ap- 
proximate the population probability by the Boltzmann 
statistical factor, ex.-p{—Ef^/T). The /3-decay rate is then 
proportional to 



g-(£.-£„)/r X Ee^E^,-ml{E, -E*^- E,)\ (22) 

which has a local maximum at Ei = 2T -\- E^ + Eq. 
That is, the /3-decay rate is large when the condition 
Ef"" = Ei - Eq 2T ^ Ee - {Eo - E^) is satisfied. On 
the other hand, it is hindered considerably in the case of 

Er<2T + ^e-(^o-^^). 

In order to investigate the temperature dependence of 
Ef^, we plot the average pairing gaps (A) in Figure [6] 
and the unperturbed energy of the first excited state 
Ef^ = Ei—Eo, evaluated with the two-quasi-particle en- 
ergy in the BCS approximation, in Figure (3 The top and 
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FIG. 5: (Color online) The /3-decay half-life T1/2 normalized to that at zero-temperature, The solid and the dashed lines 

show the results with the SLy5 and SkM* parameter sets, respectively. 



bottom panels show the results of the SLy5 and SkM* 
parameter sets, respectively. One sees that the pairing 
gaps begin to decrease significantly at temperatures of 
about T = 0.50 — 0.65 MeV {i.e., the pairing phase tran- 
sition). Likewise, the energy of the first excited state, 
also decreases rapidly at similar temperatures. For 
^^^Ru and ^^^Zr, it eventually becomes less than 2T at 
high temperatures. This should be intimately related to 
the increase of the -decay half-lives for these nuclei at 
high temperatures. On the other hand, Ef^ for ^^^Cd is 
much less sensitive to the temperature, as this nucleus is 
in the neighborhood of the double magic nucleus ^^^Sn. 



As a consequence, its /3-decay half-life monotonically de- 
creases as a function of temperature. Note that the crit- 
ical temperature for the pairing phase transition is lower 
for SkM* as compared to SLy5. This fact leads to the 
result that the /3-decay half-lives start increasing earlier 
for SkM* compared to SLy5, as can be seen in Figure [H 



IV. CONCLUSION 

We have assessed the thermal effects on /^-decay half- 
lives with astrophysical interests for even-even isotones 
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FIG. 6: (Color online) The average proton pairing gap as 
a function of temperature for the ^^°Cd, ^^^Ru, and ^^^Zr 
nuclei. The top and bottom panels are results for the SLy5 
and SkM* parameter sets, respectively. 
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FIG. 7: (Color online) Same as Fig. [6l but for the unper- 
turbed energy of the first excited state, estimated in the BCS 
approximation. 



with the neutron magic number A/" = 82. For this pur- 
pose, we have adopted the finite temperature QRPA 
method on top of the finite temperature Skyrme-Hartree- 
Fock+BCS method. We have used the to and ^3 terms 
of the Skyrme force for the particle-hole residual inter- 
action, and a J-type interaction for the proton-neutron 
particle-particle channel in the QRPA formalism. 

We have calculated the Gamow- Teller strengths in the 
temperature range from T = 0.0 to 0.8 MeV. At finite 
temperatures, new peaks appear in the strength func- 
tion due to the transitions from the excited states. From 
the calculated Gamow- Teller strengths, we have evalu- 
ated the /3-decay half-lives. As the temperature increases, 
the /3-decay half-life decreases gradually for all the nu- 
clei which we have studied. We have also found that 
the temperature dependence appears more strongly for 
nuclei with a larger atomic number. We have argued 
that this is related to the number of GT peaks in the 
strength function, determined mainly by the difference 
between the proton and neutron Fermi surfaces. We have 
also found that the /3-decay half-life begins to increase at 



T > 0.6 — 0.7 MeV for open-shell nuclei as a consequence 
of a peculiar temperature dependence of the energy of the 
first excited state due to the pairing phase transition. 

From our results, we conclude that the thermal effect 
on the /3-decay half-life is negligible at the standard r- 
process temperature, which is considered to be approxi- 
mately less than 0.2 MeV, at least for even-even = 82 
isotones. It would be an interesting future problem to 
extend the present calculations to odd-mass nuclei, in 
which the energy of the first excited state is in general 
much smaller than that in even-even nuclei and thus a 
larger thermal effects may be expected. 
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